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We study a quantum double model whose degrees of freedom are Ising anyons. The terms of the 
Hamiltonian of this system give rise to a competition between single and double topologies. By 
studying the energy spectra of the Hamiltonian at different values of the coupling constants, we 
find extended gapless regions which include a large number of critical points described by conformal 
field theories with central charge c — 1. These theories are part of the Z2 orbifold of the bosonic 
theory compactified on a circle. We observe that the Hilbert space of our anyonic model can be 
associated with extended Dynkin diagrams of affine Lie algebras which yields exact solutions at 
some critical points. In certain special regimes, our model corresponds to the Hamiltonian limit 
of the Ashkin- Teller model, and hence integrability over a wide range of coupling parameters is 
established. 
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I. INTRODUCTION 

There has been considerable interest in emergent par- 
ticles with fractional statistics, so-called anyons 1 2 . Most 
prominently, anyons appear as quasiparticle excitations 
of the ground state in the fractional quantum Hall (FQH) 
liquids 34 . Anyons are also realized in quantum spin 
models in two spatial dimensions, such as the toric 
code model 5 , the quantum dimer model on non-bipartite 
lattices^, and Kitaev's honeycomb modeP. The toric 
code model is a special case of a whole set of time- 
reversal and parity invariant lattice models that real- 
ize doubled topological quantum field theories in (2+1) 
dimension d 8 * 9 * 10 * ('quantum double models'). The im- 
plementation ofsuch models in terms of lattice spin 
Hamiltoniane P ' 13 ' 14 ' -' , or Josephson junction arrays 1 ^ is 
under active investigation. 

In this work, we are interested in the physics of a quan- 
tum double model whose microscopic degrees of freedom 
are non-abelian anyons. More specifically, we investi- 
gate a quantum double model whose degrees of freedom 
are Ising anyons. A simple example of a quantum dou- 
ble model appeared irPEl w here the degrees of freedom 
are Fibonacci anyons^ located on the links of a hon- 
eycomb lattice. The Hamiltonian penalizes Fibonacci 
anyon fluxes through the plaquettes of the lattice, and 
it is exactly solvable. This model has recently been con- 
sidered on a ladder basis, where a competing term, which 
penalizes Fibonacci anyons on the rungs of the ladder ba- 
sis, was added to the Hamiltonian , see Fig.[l] In fact, 
the system studied is an example of a quantum dou- 
ble model of non-abelian anyons with both a 'string-net 
kinetic energy (plaquette fluxes) and a 'string- net ten- 
sion' (rung fluxes), as envisioned, but not studied, irP. 
It was found that the competition between the rung and 
plaquette fluxes can be translated into the competition 
between two extreme topologies, each of them associated 
with a gapped phase. At equal magnitude of the cou- 
pling constants of rung and plaquette term, a critical 
point separating the gapped phases was observed. This 



critical point, and a second critical phase, are described 
by certain conformal field theories. The Hilbert space of 
the model is associated with a Dq Dynkin diagram which 
yields exact solutions at two critical points. 

In this paper, we focus on degrees of freedom corre- 
sponding to Ising anyons. Ising anyons are currently the 
most promising class of non-abelian anyons in the ex- 
perimental context. A fractional quasiparticle charge of 
one quarter of the electron charge (as expected for the 
Ising anyon) has been measured^ which raises hopes that 
the quasiparticles in the fractional quantum Hall state 
with filling fraction v = 5/2 are indeed Ising anyons, 
as predicted in 19 . Further systems with emergent Ising 
anyons are p + ip superconductors (or superfluids) 19 20 , 
and a quantum spin lattice modeP. The Ising theory 
possesses an additional anyon species as compared to the 
Fibonacci theory. There are two types of particle species 
(the Ising anyon a, and the fermion Hence our model 
has two coupling parameters that can be tuned: One of 
the coupling parameters varies the relative strength of 
rung and plaquette fluxes, while the other coupling pa- 
rameter varies the relative strength of the Ising anyon 
and the fermion (rung and plaquette) fluxes. We study 
the phase diagram as a function of the two coupling con- 
stants using exact diagonalization and analytical meth- 
ods. We observe extended gapped and gapless phases, 
where the latter includes a number of critical theories 
which are described by two-dimensional rational confor- 
mal field theories (rCFT) with central charges c = 1. 
These critical points are part of the Z2 orbifold of the 
bosonic theory compactified on a circle of radius R = y/2p 
(where each integer p > gives rise to a separate rCFT) . 
The Hilbert space of our model is associated with the 
extended Dynkin diagrams D4 and Dq which yields ex- 
act solutions a certain critical points. Further results are 
established by identifying certain regimes of our model 
with the quantum Ashkin- Teller model. 

The organization of this paper is as follows. After a 
brief review of some essential properties of Ising anyons in 



section II A we introduce our model in sections II B and 





c 2 




c 4 




FIG. 2: Labeling of the basis (the fusion diagram). Periodic 
boundary conditions are applied, i.e., a\ — cl2l+i, bi = &2L+1, 
and c 2 = C2L+2. 



FIG. 1: The topology associated with our model of Ising 
anyons is a high-genus surface. The ladder skeleton inside 
the sphere is a possible basis choice. In this basis, the two 
terms in the Hamiltonian project onto the 'flux' [1 (no flux), 
a, or ip] through the plaquettes (i.e., the holes of the high- 
genus surface) and the 'flux' on the rungs of the ladder basis, 
respectively, as indicated by the red arrows. 



II C We then present the results of an exact diagonaliza- 



tion study of the Hamiltonian in section III Thereafter, 
we present exact solutions of our model at certain critical 
points, based on the association of the Hilbert space with 
certain extended Dynkin diagrams (section IV). 



II. MODEL 
A. Ising anyons 

In the following, we recapitulate some essential prop- 
erties of the degrees of freedom of our model, so-called 
Ising anyons^. There are three different particle 'species' 
in the Ising theory, the trivial particle 1, the Ising anyon 
(j, and the fermion ip. 

The coupling of two Ising anyons is determined by the 
fusion rules, which are the analogs of Clebsch- Gordon 
rules for ordinary angular momenta. The fusion rules of 
the Ising theory are given by 



a fourth anyon d, 
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Here, labels a, 6,..., take values 1, cr, and ip, and the di- 
agrams represent the quantum states of the 'four anyon 
system' where anyons fuse in the specified order. The 
non-trivial elements (F£ bc )f (i.e., (F^ bc )( / 1) of the 
Ising theory are (F^ a )° = -1, and 
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The modular S'-matrix is a basis transformation which 
relates the anyon 'flux' of species b through an anyon loop 
of species a to the case without anyon loop by 
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for the case of Ising anyons. 



C rXcr = l+'0, crX^ = a, ^X^ = l, 1x1 = 1. (1) 

These fusion rules can be written in terms of the fusion 
matrices Nj whose entries (Nj) 3 ^ equal to one iff the 
fusion of anyons of types ji and 22 into j is possible. The 
fusion rules are related to the quantum dimensions dj, 
j = 1,(7,^, by Njd = djd, where d is the eigenvector 
corresponding to the largest positive eigenvalue of the 
matrix Nj. The quantum dimensions of the Ising theory 
are di = 1, d a = y/2 and d^ = 1, and the total quantum 



dimension is V = 




In analogy to the 6j-symbols for ordinary SU(2) spins, 
there exists a basis transformation F that relates the 
two differents ways three anyons a, 6, and c can fuse to 



B. Hilbert space 

Anyonic degrees of freedom are non-local, i.e., the 
Hilbert space of a multi-anyon system is not the ten- 
sor product space of Hilbert spaces associated with lo- 
cal degrees of freedom, as is the case for ordinary spins. 
The Hilbert space of a multi-anyon system can be rep- 
resented in terms of a fusion diagram which is a triva- 
lent graph with each line segment symbolizing a certain 
anyon species, and the fusion rules being obeyed at the 
vertices (see Fig. [2|. Each distinct occupation of the 
fusion diagram represents a basis state, and the inner 
product of two identical states is one, while the inner 
product of two different states is zero. By means of 
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FIG. 3: A different decomposition of the high- genus surface 
into three-puncture spheres (as compared to the one shown 
in Fig.[T]), yields a different choice of basis (the black graph), 
also shown in Fig. [2] 



(for example) F-transformations, different basis choices 
of the same system can be related. It is the topology 
(here, the high-genus surface of Figs. [I] and [3| that defines 
the Hilbert space. Different basis choices correspond to 
different decompositions of the high-genus surface into 
three-punctured spheres, as can be seen by comparing 
Figs, [l] and [3] We formulate our Hamiltonian in the ba- 
sis choice which is shown in Figs. [2] and [3] 

In the terminology of Fig. [2j the occupations of ele- 
ments Ci are either q = a (i = 2, 4, 2L), or q G {1, ^} 
(i = 2,4,..., 2L). This means that there are two indepen- 
dent sectors of the Hilbert space of our model: 

• IS (Integer sector): a G {1,^} (« = 2,4, ...,2L), 
(a u h) G {(l,l),(a,a),(^^),(^l),(l,^)} (i = 
1,3,...2L- 1). 

• HIS (Half-integer sector): q = a (i = 

2,4,..., 2L), {(LiM) £ {(1,(7), ((7,1), (^,(7), ((7,^)} 

(i = 1,3,...,2L- 1). 

Using the fusion matrices Nj, it is straightforward to 
evaluate the number of basis states, 5, as a function of 
the number of plaquettes, L. We apply periodic bound- 
ary conditions, i.e., a x = a 2 L+i, h = &2L+1, c 2 = c 2 l+2- 
The number of basis states is given by 

B = £ (^^(^^(^^-(^J^ 1 (6) 
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where the summation | a . 6 . c i runs over all possible 
labelings of the basis. 
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FIG. 4: F-transformation of a local element of the basis of 
Fig. [2] to the ladder basis. 



C. Hamiltonian 

The Hamiltonian contains two non-commuting terms 
which act in alternating manner on even and odd labels 
z, i = 1,2,...,2L (terminology as in Fig. [2). The pla- 

quette operator P^ projects onto anyon flux s (where 
s G {1, cr, iJj}) through a plaquette indexed by an odd in- 
teger i. The rung operator r\ s ^ (i even) projects onto 
the anyon occupation of a rung (i.e., it is diagonal in the 
ladder basis of Fig. [I] but not in the basis used [Fig. [2]). 
In the most general form, the Hamiltonian is given by 
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We set J p = cos(0), J r = sin(0), J\ = cos(0), 
sin(0) and rewrite the Hamiltonian as (note that P- 
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The parameter # controls the dimerization of the model. 
If J r = J p , i.e., 9 = 7r/4 or = 57r/4, the dimerization 
is zero, i.e., the local terms Hi (where Hi = Pi if i odd, 
and Hi = Ri if i even) have identical coupling strengths 
at each 'site' i. 

The exact form of the terms P^ and R^ was dis- 
cussed irff^ for the case of Fibonacci anyon degrees of 
freedom. However, we shall repeat the derivation of this 
non-standard Hamiltonian for the case of Ising anyons. 

(s) 

We begin with the local plaquette term P> \ We insert 
an additional anyon loop of type t G {l,cr, ^} into the 
center of the plaquette composed by variables (a^,^), 
and project onto the flux through this additional loop 
(and hence the flux through the plaquette) using the S- 
matrix Eq. Q, 
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We proceed further as follows, 




^1 (^Oiait)l' (^c/ + i6it) a i 



(-^ai ait ) 1* (^Ci+i 6^ ) a i (^Ci*_ i a* t ) 

a/ ,b'.,m 




~ ( ^Ci + 1 b % t ) a i ( ^cL i a 4 1 ) b\ 

where we used the identity 
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and the orthogonality relation T, e ( F abc) e f( F dab)e = $e,k- 
By using an F-transformation, it is possible to trans- 
form between the basis of Fig. [2] and the ladder basis, as 
shown for a local element in Fig. [4j Using such a trans- 
formation, the projector onto a rung with occupation s 
is given by 

Ri ^\ C i) ~ 5Z(^i-ia*-ia i+ i)ci (^6i^i ai_ i a i+i )c / i I C i) ■ 

(11) 

It is straightforward to construct a matrix represen- 
tation of the Hamiltonian (|7|). In the half-integer sec- 
tor (HIS), the variables at even sites i are fixed, i.e., 
Ci = 0. We associate the local 'site' variables (a^b^) £ 
{(1,(7), ((7,1), (^cr), (o", VO} (* °dd) with the four 
unit vectors in four dimensions, respectively, and de- 
fine = Diag(l, 0,0,0), = Diag(0, 1, 0, 0), 
n W,<r) = Diag(0, 0,1,0), and = Diag(0, 0, 0, 1). 
Evaluating Eqs. ([9| a nd (pT| using the F- and ^-matrix 
elements (see section |IIA| ) yields a 4 x 4 representation 
of the Hamiltonian in the HI sector, 
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where 
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The Hamiltonian (12) is invariant under variable exchanges (1,0) <-> (^,(7) and (0, 1) <-> (0, ^) (independently and 
simultaneously). It is also invariant under simultaneous exchange of variables (1,0") <-> (0, 1) and (ip,(r) <-> (0, -0). 

In a similar manner as for the half-integer sector, it is possible to construct a 7x 7 matrix representation of the Hamil- 
tonian ([7]) in the integer sector. Let the variables 1, (1, 1), (1, -0), (0, 0), (ip, 1), (-0, -0), -0 be associated with the seven 

^ ~ m^ 1 ' 1 ) = Diag(0,l,0,0,0,0,0), 



unit vectors in seven dimensions (in this order) and let m 1 = Diag(l, 0, 0, 0, 0, 0, 0), 
and so on. In this notation, the Hamiltonian Q in the integer sector takes the form 
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Hamiltonian (13) is invariant under variables exchanges 
(i) (1,1) <-> (0^0) an d ( n ) (l?^) ^ (^,1) (indepen- 
dently and simultaneously). It is also invariant under 
simultaneous exchange of 1 <-> 0, (^,1) <-> (0,0) and 
(l,^) <-> (1,1) (where exchanges (i) and (ii) are applica- 
ble, too). 



D. The half-integer sector and the quantum 
Ashkin- Teller model 



III. NUMERICAL RESULTS 

We first outline the topological feature of our model 
that determines its criticality at equal magnitude of rung 
and plaquette term. We recapitulate the identification of 
a conformal field theory based on the energy spectrum in 
a system of finite size and review the operator content of 
the Z2 orbifold of the compactified bosonic theory. Then, 
we present the results of the exact diagonalization of the 
Hamiltonian matrix. 



The half-integer sector of the high-genus ladder of Ising 
anyons studied in this paper is equivalent to the quantum 
Ashkin- Teller model. The quantum Ashkin- Teller model, 
which can also be mapped onto the staggered XXZ chain, 
is the Hamiltonian limit of the classical Ashkin- Teller 
model, and it was first studied in 22 . 

The correspondence of our model in the half-integer 
sector and the quantum Ashkin- Teller model becomes 
immediately apparent when relabeling the indices of 
Hamiltonian Eq. (12) according to i—1 — > n, z+1 — > n+1, 



and comparing this Hamiltonian with the Hamiltonian 
given irPH (see alscP^). Relating the coupling constants 
6 and to the ones irP^ allows us to confirm the results 
of sections HID an dJIV A , and add further details to the 
phase diagram Fig. [5] We believe that the numerical re- 
sults discussed section |IIID| are of interest despite prior 
studies of the quantum Ashkin- Teller model, and we note 
that the derivation of the exact solution in section HV Al 
is a consequence of the unique structure of the Hilbert 
space of our anyonic model. 



E. Numerical method 

We diagonalize the Hamiltonian matrix using the 
Lanczos algorithm 24 . By employing periodic boundary 
conditions, we obtain the energy eigenvalues as a func- 
tion of momenta k x = 27rn/L, n = 1,2, as well as 
k y = 0, 7r (invariance of the Hamiltonian under exchange 
of the di and bi variables, this symmetry corresponds 
to simultaneous variable exchanges (l,cr) <-> (a, 1) and 
(0, a) <-> (a, in the half-integer sector, and exchange 
(l,^) <-> (0>, 1) m the integer sector). We employ an 
implementation of the Lanczos algorithm in the ALPS 
librar}E^. 



A. Competing topologies 

The competition between the rung and plaquette terms 
correspond to a competition between single and double 
topologies 17 . This can be understood by switching to the 
ladder basis of Fig. [I] In the ladder basis, the plaquette 
term projects onto the flux s through the plaquette, and 
the rung term projects onto the flux s on the rung of 
the ladder. We consider the Hamiltonian at the points 
J p = 1, J r = 0, J-0 = {9 = 0, (j) = 0), and J r = 1, 
J p = 0, J0 =0 (6 = 7r/2, = 0), respectively. For the 
former choice of coupling constants, the rung term is zero, 
and the Hamiltonian favors the absence of a- and 0-fluxes 
through the plaquettes. However, if there are no fluxes 
through the holes of the high- genus surface, the holes can 
be closed, and we are left with a single cylinder (a torus 
for the case of peridic boundary conditions). In contrast, 
at the latter choice of coupling parameters, the plaquette 
term is zero, the Hamiltonian favors the absence of a and 
particles on the rungs. Hence, the rungs can be 'cut 
off', and the resulting surface is that of two independent 
cylinders (two tori for periodic boundaries). In this work, 
we mainly consider the points of equal magnitude of rung 
and plaquette term where the competition between single 
and double topologies renders the system critical over a 
large range of coupling parameters. 



B. Identification of conformal field theories 

The spectrum of a conformal field theory (CFT) in a 
system of finite size L and periodic boundary conditions 
has the following energy eigenvalue^ 
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TABLE I: Scaling dimensions /io + (aside from the ground 
state with h — h = 0) of the operators of the Z2 orbifold of the 
boson compactified on a circle of radius R = \/2p for some 

theories (i.e., some integer p > 0). The scaling dimensions 

- 2 

which depend on p are given by h + /io = where n = 
1, 1, and by p/2. The following abbreviations are used: 

scft = superconformal CFT with c = 1, Potts = 4-state Potts 
theory, para = parafermion CFT with c = 1, Ising 2 = square 
of the Ising CFT, KT = Kosterlitz-Thouless transition. 



where c is the central charge of the CFT, and the velocity 
v is an overall scale factor. The scaling dimensions h + h 
take the form h = h° + n, h = h° + n, where n and n 
are non- negative integers, and h° and /z° are the holo- 
morphic and antiholomorphic conformal weights of pri- 
mary fields of a given CFT of central charge c. Energies 
with h and h such that n and ft zero are associated with 
primary fields while energies with n and/or n non-zero 
correspond to descendant fields. There are some con- 
straints on the momenta k x (in units 2ir / L): k x = h — h 
or k x = h — h + L/2. The system size L corresponds to 
the number of plaquettes of the basis, also denoted by 
L in the previous and upcoming discussion. By rescaling 
the eigenenergies obtained from exact diagonalization ac- 
cording to Eq. (14), we are able to identify a number of 



conformal field theories. 



C. Z2 orbifold of the boson compactified on a circle 
of radius R = y/2p 

As was mentioned in the introduction, we identify a 
number of conformal field theories with central charges 



1CDP 



FIG. 5: Half-integer sector (HIS): Phase diagram at angles 
6 = 7r/4, and <fi G [0, 2tt) 32 . The coupling constants are Ji = 
cos(0) and = sin (</>). The positions of some of the c = 1 
theories (boson orbifold compatified on a circle of radius R = 
y/2p) are indicated, see table [i] for details on the operator 
content. For example, p = 4 stands for the 4-state Potts 
CFT. 



c = 1. We observe that these theories are part of the Z2 
orbifold the the bosonic theory compactified on a circle 
of radius R = ^/2p, w ith e ach integer parameter p > 1 
defining a rational CFT^EI] Aside from the ground state 
(h = h = 0), there are two fields with scaling dimension 
ho + ho = 1/8, two fields with scaling dimension 9/8 
(these four operators are the so-called twist operators), 
one field with scaling dimension 2, two fields with scaling 
dimension p/2, and p — 1 fields with scaling dimensions 
n 2 /2p, n = 1,2, — 1 (see Table[l|. The more promi- 
nent of the critical theories of the Z2 boson orbifold are 
the Kosterlitz-Thouless theory (p = 1), the theory of two 
decoupled Ising models (p = 2), the c = 1 parafermion 
CFT (p = 3), the 4-state Potts model (p = 4), and the 
c = 1 superconformal CFT (p = 6). A number of the 
orbifold theories (those with p < 4) are observed in a 
critical line of the Ashkin- Teller model 29 which is a two- 
dimensional classical lattice model of two decoupled Ising 
models which are coupled by a four-spin interaction. 

There exists a relation between the c = 1 orbifold 
theories and the extended Dynkin diagrams D n of the 
simply-laced affine Lie-algebras of type D: for p = m 2 , 
m = 1, 2, the corresponding Z2 orbifold theory is asso- 
ciated with the extended Dynkin diagram '. The 

extended Dynkin diagrams D n define so-called restricted- 
solid-on-solid (RSOS) models which are 2D statistical lat- 
tice models whose degrees of freedom are integer-valued 
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FIG. 6: HIS: Rescaled energy spectrum (from exact diagonal- 
ization at system size L — 10) at = 7r/4, cj) = 7r/2, and the 
CFT assignments of the 4-state Potts model. There are three 
fields with scaling dimension 2 (see Table 0, at momentum 
k x = 0, however, the finite-size effects are rather strong. 



FIG. 7: HIS: Rescaled energy spectrum (from exact diagonal- 
ization) at = 7r/4, cj) = 0.557T, and the CFT assignments of 
the parafermion CFT with c — 1. 



heights on the nodes of the lattice with the constraint 
that heights on nearest-neighbouring lattice sites are ad- 
jacent nodes in the defining Dynkin diagram. The parti- 
tion function of these RSOS models is a discrete version 
of the partition function of the rCFTs associated with 
the respective Dynkin diagrani 3Q * 31 i 



D. Numerical results in the half-integer sector 



(HIS): \J r \ = \J P 



In this section we discuss the results of the exact di- 



agonalization of the Hamiltonian (12) for equal magni- 
tude of plaquette and rung terms, i.e., \J p \ = \J r \. In 
the following, we refer to the case = 7r/4, however, all 
other cases of equal strength of rung and plaquette cou- 
pling (0 = 3tt/4, = 5tt/4, = 7tt/4) yield the identical 
results^. 



1. Gapless phases 

At angles = and = 7r/2, the model is critical and 
described by the 4-state Potts model, as shown in Figj6j 
which is confirmed by an exact solution (see section |Tv| ). 
In fact, for all angles G [— tt/2,0] and G [tt/2, 7r]the 
system is gapless with continuously varying critical ex- 
ponents. We are able to match the energy spectra at dif- 
ferent angles to several of the orbifold theories, as indi- 
cated in the phase diagram Fig. [5] (see also a figure of the 
c = 1 parafermion CFT Fig. [7]) . The two gapless phases 
in the phase diagram (Fig. [5]) differ by the /^-momentum 



quantum numbers of the twist operators. In one of 
the gapless phases, the four twist operators with scal- 
ing dimensions 1/8, 1/8, 9/8, 9/8 have momentum quan- 
tum numbers (k xi k y ) = (0, 0), (0, 7r), (0, 0), (0, 7r), while 
in the other gapless phase the momenta are (k x ,k y ) = 
(-7T, 0), (tt, 7r), (7t, 0), (7r, 7t). The eigenenergies associated 
with the remaining operators always appear in momen- 
tum sector k x = 0. The fields with scaling dimensions 
n 2 /2p, n = 1,. — 1, have momentum quantum num- 
bers k y = if n is even, and k y = 7r if n is odd. The 
numerical results indicate that fields with scaling dimen- 
sions p/2 are both in momentum sector k y = for p 
even, while for p odd, these two fields have momentum 
numbers k y = and k y = n, respectively. The marginal 
operator is in momentum sector (k x ,k y ) = (0,0). We 
located the critical theories at angles = and <j) = 0.57T 
(p = 4), (j) w 0.55tt and w 1.95tt (p = 3), (j) w 0.75tt 
and (j) w 1.757T (p = 2), and <p « 0.957T and <p « 1.557T 
(p = 1). The exact locations of these orbifold theories 
can be determined by comparison wittP. The Hamil- 
tonian (12) is integrable along the orbifold line, i.e., be- 



tween the Kosterlitz-Thouless point and the 4-state Potts 
point (see^and references therein). 

The region between the Kosterlitz-Thouless transition 
(p = 1) and the points (j) = tt and (j) = 37r/2, respec- 
tively, is gapless with continuously varying critical ex- 
ponents. This region is denoted as 'critical fan' in 22 . 
The transition points <p = tt and <p = 37r/2 are gap- 
less, and exhibit a three- fold degenerate ground state 
[momenta (k Xl k y ) — (0, 0), (jr, 0), (jr, 7r) at <j) = 7r, and 
(k x ,k y ) = (0,0),(0,7r),(7r,7r) at = 3tt/2]. They mark 
a first order transition between the critical phases and 
gapped phase I, as indicated by a jump of the ground 
state energies (as a function of (j)) at these angles (not 
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shown here). These two critical points are in the uni- 
versality class of the potassium dihydrogen phosphate 
(KDP) modeP. 



2. Gapped phases 

We briefly discuss the two gapped phases in the phase 
diagram Fig. [5] In gapped phase I [0 G (7r,37r/2)], the 
ground state is two-fold degenerate. Above the ground 
state, a flat quasiparticle band is observed. At angle 
= 57r/4, the Hamiltonian is of form H = — | P^li — 

\ ^2 i R^2i ) • At this point, the energy is minimized if all 
rungs have occupation cr. This is realized for any config- 
uration of form [we omit the indices q = cr, i.e., = 
|(ai,6i),(a 3 ,6 3 ),...)] = \(a u cr), (cr, 63), (a 5 , a), ...), 

where ai, 63,(25,... £ {1,^}. All states of this form, 
and hence also the ground states, appear only in mo- 
mentum sectors (k xi k y ) = (0,0) and (k xi k y ) = (7r,7r). 
The ground states are the product states of local states of 
form \{a % A)) = ^(|(1, <x)> - *))) (i = 1, 5, ...2L-3), 
and 1(^,60) = ^(l(M)H(<^)>) (* = 3,7,...,2L-1), 
and they are hence a superposition of all states of form 
l^j), where the magnitude of the weights depends on the 
multiplicities of the states according to the symmetries. 
The numerical results confirm that this is indeed the cor- 
rect construction for any point in the gapped phase I. 
The symmetry and the exact form of the ground states 
is indicative of a system of two independent sublattices. 
Gapped phase I corresponds to the 'antiferromagnetic 
frozen phase' irPH 

In gapped phase II, the ground state is also two-fold 
degenerate, and the quasiparticle dispersion has a lead- 
ing cosine shape. At coupling parameter = 7r/4, the 



At 



Hamiltonian is of form H = | P 2 i-i + \ R' 
this point, the energy is minimized if all rungs (in the lad- 
der basis) have occupation 1 or ip. This is realized for any 
configuration of form \^n) = |(<2i, cr), (03, cr), (as, cr), ...) 
where ai, (23,(25, ... G {l,^}. All states of this form, and 
hence also the ground states, appear only in momen- 
tum sectors (k Xl k y ) = (0,0) and (k Xl k y ) — (0,7r). The 
ground states are the product states of local states of 
form 1(^,6*)) = ^(|(l,a)) + |(^cr))), i = 1,3,...,2L-1 
and thus are a superposition of all states of form \^n) 
where the magnitude of the weights depends on the mul- 
tiplicities of the states according to the symmetries. This 
gapped phase corresponds to a partially ordered phase 
(ordered in one of the two Ising spins) 



E. Numerical results in the integer sector (IS): 



Ur 



\Jv 



We discuss the results of the exact diagonalization of 
the Hamiltonian ( fl3| ) for equal magnitude of plaquette 
and rung term. In the following we refer to the case 

= n/m 




ex. sol, 
E=0,0.5,l,1.5 ; 



FIG. 8: Integer sector (IS): Phase diagram at angles = 7r/4, 
and G [0, 27rj^. The coupling constants are Ji = cos(0) 
and — sin(0). The positions of some of the c — 1 theories 
(boson orbifold compatifled on a circle of radius R — y/2p) 
are indicated, see table [I] for details on the operator content. 
There exist exact solutions (ex. sol.) at seve ral po ints: = 0, 
4> = tv/2 (p = 16 boson orbifold, section IV B), = 7r/4, 
cj) = 5tt/4 (section |IIIE2 ) . 



1. Gapless phases 

The phase diagram in the integer sector (IS) at equal 
strength of rung and plaquette term is similar to the 
one of the half-integer sector, with two extended gapped 
and two extended gapless phases, as illustrated in Fig. [8] 
However, a much larger number of the rational c = 1 
critical theories appears. At angles = and <p = tt/2 
(where we found the p = 4 theory in the half-integer 
sector), we identify the p = 16 orbifold theory, as can 
be seen in Fig. [9] This numerical result is confirmed by 
the exact solution given in section IV B We identify the 
c = 1 orbifold theories for p = 1, 2, 16 in the two gap- 
less phases. One example of the observed theories is the 
c = 1 superconformal CFT (p = 6) which is shown in 
Fig. [lOj As was the case for the half-integer sector, the 
two gapless phases differ by the k x momentum quantum 
numbers of the twist operators. The momentum assign- 
ments of the twist operators in the integer sector are 
identical to the ones in the half-integer sector. In con- 
trast to the half-integer sector, all remaining fields have 
momentum quantum numbers (k Xl k y ) = (0,0). 

The angles at which some of the critical theories are 
located are the following: = and = 0.57T (p = 16) 
p = 16, w 0.857T and ~ 1.657T (p = 6), ~ 0.957T and 
w 1.55tt (p = 4), w 0.97tt and w 1.53tt (p = 3), 
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FIG. 9: IS: Rescaled energy spectrum (from exact diagonal- 
ization) at = 7r/4, cj) = tt/2, and the CFT assignments of 
the Z2 orbifold of the compactified bosonic theory at radius 
R ^/2p with p 16. 



FIG. 10: IS: Rescaled energy spectrum (from exact diagonal- 
ization) at = 7r/4, cj) = 0.857T, and the CFT assignments of 
the Z2 orbifold of the compactified bosonic theory at radius 
R — ^J2p with p = 6 (c = 1 superconformal CFT). 



<j> « 0.985tt and <j> w 1.515tt (p = 2), and « 0.995tt and 
0^ 1.505tt (p= 1). 

It is difficult to determine the exact position of the 
transition between either of the gapless phases and the 
gapped phase I (dashed lines in Fig.|8|. However, the en- 
ergy eigenvalue associated with the field with scaling di- 
mension 9/2p does not become degenerate with the eigen- 
value associated with the twist fields of scaling dimension 
1/8 when approaching the gapped phase I (from either 
side). This means that the orbifold theory with p = 36 
does not appear, and thus the orbifold theory with p 
must be one of the theories with 16 < p < 36. From 
the analyti cal re sults applicable to the half-integer sec- 
tor (section II D ) it is known that the transition between 
gapped phase I and the critical phases is located at <p = 
and (j) = 7r/2, respectively. It therefore seems likely that 
the corresponding transition in the integer sector is also 
located at these angles. 



2. Gapped phases 

We very briefly remark on the gapped phases of phase 
diagram Fig. [8j In gapped phase I, a flat quasiparticle 
band is observed above a highly degenerate ground state. 
These degenerate ground states are superpositions of 
states where local variables (a^bi) = (cr, a) are followed 
by variables (a i+2 ,6 i+2 ) = (1,1), = 
K+2A+2) = 0,1) or (a i+2 ,6 i+2 ) = (1,^), and vice 
versa. In gapped phase II, the ground state is non- 
degenerate, and the quasiparticle band exhibits a leading 
cosine shape. 



F. Phase diagram for \J r \ 7^ \J P 



The equivalence of the half-integer sector model (12) 



with the Ashkin- Teller model yields the phase diagram of 
the half-integer sector (for details of the various phases 
in spin language see^). We numerically verified the ex- 
istence of a critical point which is described by the Ising 
CFT with central charge c = 1/2 for coupling param- 
eters = 0.32-7T and <p = 0.3527T, which is in agree- 
ment with prior results on the quantum Ashkin-Teller 
modeP^^. Each of the three primary fields of the Ising 
theory appears twice, where momentum symmetry sec- 
tors are (k x , k y ) = (0, 0), (0, n) for operators with scaling 
dimensions and 1, and (k x , k y ) — (7r, 0), (7r, tt) for oper- 
ators with scaling dimension 1/8. 

In the integer sector, the exact diagonalization of 
Hamiltonian (13) at coupling parameters 6 = 0.327T, 



cj) = 0.3527T yields that the corresponding phase is of 
Ising universality, too. However, the field with scaling 
dimension 1/8 appears with triple degeneracy [momenta 
(k x , k y = (0, 0), (7r, 0), (7r, 7r)], while the fields with scaling 
dimensions and 1 appear only once [both at momentum 
(k x ,k y ) = (0,0)]. Further details of the phase diagram 
in the integer sector away from \J r \ 
studied. 



J p \ remain to be 



IV. EXACT SOLUTIONS AT CRITICAL 
POINTS: DYNKIN DIAGRAMS 

In this section, we identify the Hilbert space of our 
model with certain extended Dynkin diagrams. We ob- 
serve that the Hamiltonian, at certain coupling con- 
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FIG. 11: D4 Dynkin diagram associated with the Hilbert 
space of the h alf integer sector (HIS). The symmetries of the 
Hamiltonian (12) under exchange of variables are indicated 
by green arrows. 



WM) 



FIG. 12: Dq Dynkin diagram associated with the Hilbert 
space o f the integer sector (IS). The symmetries of the Hamil- 
tonian ( 13 ) under exchange of labels are indicated by green 
arrows. 




stants, corresponds to the restricted-solid-on-solid mod- 
els which are associated with these extended Dynkin di- 
agrams. 



A. Half-integer sector (HIS) 

We associate a label q = cr (terminology of Fig.[2| with 
the even numbered 'sites' i. With the odd- numbered 
'sites' i we associate a variable consisting of a pair 
of labels, (a^,^), which can assume four values, i.e., 
(a,i,bi) = (l,cr), (a,i,bi) = ((7,1), (a,i,bi) = and 
(ai,bi) = (cr, t/>). If variables (a^b^) and ci±\ are allowed 
to meet at the vertices (as a consequence of the fusion 
rules) they are adjacent nodes on the Dynkin diagram 



form a representation of the Temperley-Lieb algebra^, 



of the extended D4 Lie algebra, as illustrated in Fig. 11 
Any local label (a^,^) at an odd- numbered site i allows 
for label q_i = a at the neighboring even- numbered 
sites, which is reflected in the fact that label a is con- 
nected by a line all four possible labels [a^bij in the 



Dynkin diagram Fig. 11 



The adjancency matrix 3 -^ of the Dynkin diagram of 
Fig. [TT] is given by 



x d 4 



/ 1 \ 
10 
110 11 
10 

\ 1 ) 



(15) 



where the matrix indices are associated with the 
five different variables in the following order: 
(cr, -0), (ip, cr), cr, (cr, 1), (1, cr). The largest eigenvalue 
of is 2, and the corresponding eigenvector is given 
by 

V = {V(a^)^(xl;,a)^a,V(a,l)^(l,a)) = (1,1,2,1,1). (16) 

The operators 

ei\xi, Xi-i,Xi,Xi+i, ...x 2 l) (17) 
i( e i)xit\]xi\ x i, •.;Xi-i,X f i ,X i+1 , ...,a?2L>, 



el 



Dei , 

[ei, 6j] = for \i 



(18) 



•il >2, 



where V = 2. At coupling constants 6 = 7r/4, <p = Cp^, 
the local terms of our Hamiltonian Eq. ([7]) (in the half- 
integer sector) Hi = (i odd) and Hi — (i even) 
equal to Hi = — \^% which can be seen by evaluating 
the operators e^. It can be shown that the Hamiltonian 
(at these coupling constants) defines the two-row transfer 
matrix of the RSOS model that is associated with the 
£>4 Dynkin diagram 3 -^. Consequently, our model in the 
HI sector at angle s = 7r/4, (j) = (P2J is described by the 
4-state Potts CFT^EQl. 



B. Integer sector (IS) 



In analogy to the discussion of section [TV A[ we asso- 
ciate a label q = 1 or q = ip with the even numbered 
sites, while the odd- numbered sites are associated with 
variable consisting of a pair of labels, (a^,^) which can 
assume five values, i.e., (a^bi) = (1,1), (cbi,bi) = (cr, a), 

(diM) = (^VO, ( a iM) = W>,1) and (diM) = (1,V0- 
Variables {a^bij and Ci±i that may fuse at the vertices 
[according to the fusion rules Eq. 0] are adjacent nodes 
on the Dynkin diagram of the extended Dq Lie alge- 
bra, as illustrated in Fig. 12 For example, a local la- 
bel (ai,bi) = (cr, cr) at an odd-numbered site i allows for 
labels Ci±i = 1 and Ci±\ — ip at the neighboring even- 
numbered sites, which is reflected in the fact that label 
(cr, cr) is connected by a line to both labels 1 and ip in the 
Dynkin diagram. The components v x . of the eigenvec- 
tor associated with the largest eigenvalue of the adjan- 
cency matrix of the Dq diagram define a representation 
Eq. (17) of the Temperley-Lieb algebra associated with 
the Dq diagram. Again, it is straightforward to verify 
that the Hamiltonian in the integer sector at parameters 
= tt/4, (j> = C 32 is of form H = -\ ^ e*. This means 
that the Hamiltonian (at these coupling parameters) is 
that of the RSOS model defined by the Dq Dynkin di- 
agram, and hence the critical theory is the p = 16 Z2 
boson orbifold theor^E^. 



11 



V. CONCLUSIONS 

We study a quantum double model whose degrees of 
freedom are Ising anyons, and whose Hamiltonian imple- 
ments a competition between single and double topolo- 
gies. We observe a series of quantum critical points de- 
scribed by conformal field theories with central charge 
c = 1. These critical theories are part of the Z2 orb- 
ifold of the bosonic theory compactified on a circle. By 
associating the Hilbert space of our model with certain 
extended Dynkin diagrams, we find exact solutions of our 
model at some critical points. In one of its Hilbert space 



sectors, our model corresponds to the quantum Ashkin- 
Teller model. 

This work demonstrates the exciting physics of quan- 
tum double models which are of great interest in the con- 
text of topologically ordered phases of matter and topo- 
logical quantum computation. It contributes further to 
the unde rstanding of models of interacting non-abelian 
anyon d 17 ' 37 ' 38 ' ™ . 

The author thanks the anonymous referee for useful 
comments and E. Ardonne, D. Huse, A. Kitaev, A. Lud- 
wig, S. Trebst, M. Troyer and Z. Wang for enjoyable col- 
laboration on related work. 



1 J. M. Leinaas and J. Myrheim, // Nuovo Cimenti 37, 1 
(1977). 

2 J. Frohlich and F. Gabbiani, Rev. Math. Phys. 2, 251 

(1990) . 

3 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

4 A review of the physics of (non-abelian) anyons and the rel- 
evant physical systems can be found in: C. Nayak, et. al, 
Rev. Mod. Phys. 80, 1083 (2008). 

5 A. Yu. Kitaev, Ann. Phys. 303, 2 (2003). 

6 R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2001). 

7 A. Yu. Kitaev, Ann. Phys. 321, 2 (2006). 

8 M. Levin and X.-G. Wen, Phys. Rev. B 71, 045110 (2005). 

9 P. Fendley and E. Fradkin, Phys. Rev. B 72, 024412 (2005). 

10 M. H. Freedman, et. al, Ann. Phys. 310, 428 (2004). 

11 G. K. Brennen, M. Aguado and J. I. Cirac, New J. Phys. 
11, 053009 (2009). 

12 R. Konig, arXiv:090 L1333| [quant-ph]. 

13 F. A. Bais and J. C. Romers, Ann. Phys. 324, 1168 (2009). 

14 L. Jiang, et. al, Nature Physics 4, 482-488 (2008). 

15 B. Doucot, L. B. Ioffe and J. Vidal, Phys. Rev. B 69, 
214501 (2001). B. Doucot, et. al, Phys. Rev. B 71, 024505 
(2005). B. Doucot and L. B. Ioffe, New J. Phys. 7, 187 
(2005). L. B. Ioffe, et. al, Nature 415, 503-506 (2002). 

16 The Fibonacci theory is a simple class of non-abelian 
anyons, which contains only one anyon species, t he F i- 
bonacci anyon r. For more details, see, for example) 4 * 40 !. 

17 C. Gils, et. al, to be published in Nature Physics, 
|arXiv:0906.1579[ cond-mat]. 

18 M. Dolev, et. al, Nature Physics 452, 829-835 (2008). I. P. 
Radu, et. al, Science 320, 899-902 (2008). 

19 G. Moore, and N. Read, Nucl. Phys. B 360, 362-396 

(1991) . 

20 N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

21 For a review of the theory of Ising anyons see, for example, 
Refs™ 

22 M. Kohmoto, M. den Nijs and L. P. Kadanoff, Phys. Rev. 
5 24,5229 (1981). 

23 G. von Gehlen and V. Rittenberg, J. Phys. A 20, 227 
(1987). 

24 C. Lanczos, J. Res. Nat. Bur. Stand. 45, 255-282 (1950). 

25 http: / / alps.com-phys.org 

26 J. L. Cardy, J. Phys. A 17, L385 (1984); Nucl. Phys. B 
270, 186 (1986). 

27 P. Ginsparg, Nucl. Phys. B 295, 153 (1988). 

28 R. Dijkgraaf, C. Vafa, E. Verlinde, H. Verlinde, Commun. 



Math. Phys. 123, 485 (1989). 

29 J. Ashkin and E. Teller, Phys. Rev. 64, 178 (1943). See also 
Baxter, R. J. Exactly solved models in statistical mechan- 
ics. Academic Press, London, 1982, and references therein. 

30 V. Pasquier, Nucl. Phys. B 285, 162 (1987); J. Phys. A 
20, L221 (1987); J. Phys. A 20, L1229 (1987). 

31 A. Kuniba and T. Yajima, J. Stat. Phys. 52, 829 (1987). 

32 The choice — tt/4 and <j> G [0, 2tt) is equivalent to the 
choices = 3tt/4 and ^ 7r — <j), — bir /4 and cj) <— tt + <j), 
as well as = 7tv/4 and 4> <— 2ty — (j). 

33 E. H. Lieb, Phys. Rev. Lett. 18, 692, 1046 (1967); Phys. 
Rev. Lett. 19, 108 (1967). 

34 The adjacency matrix A is the matrix whose only non-zero 
matrix elements are A ajOC > = 1 when a and a are nearest 
neighbors on the Dynkin diagram. 

35 N. Temperley and E. Lieb, Proc. Roy. Soc. Lond. A 322, 
251 (1971). 

36 We do not present this calculation here. However, it pro- 
ceeds analogously to the one given irl 3 - 7 4 

37 A. Feiguin,et al. Phys. Rev. Lett. 98, 160409 (2007). 

38 S. Trebst, et. al, Phys. Rev. Lett. 101, 050401 (2008). 

39 C. Gils, et. al, to be published in Phys. Rev. Lett., 
arXiv:0810.2277 cond-mat]. 

40 S. Trebst, M. Troyer, Z. Wang, A.W.W. Ludwig, Prog. 
Theor. Phys. Supp. 176, 384 (2008). 



